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ABSTRACT 

'Water In Star-forming regions with HerscheV (WISH) is a key programme dedicated to studying the role of water and related species during 
the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the 
Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 
in several Hj^O, Hj^O, and CO transitions. Line profiles are resolved for five H^^O transitions in each source, revealing them to be surprisingly 
complex. The line profiles are decomposed into broad (>20 kms"'), medium-broad (~5-10 kms"'), and narrow (<5 kms"') components. The 
Hj^O emission is only detected in broad 1 lo-loi lines (>20 km s"'), indicating that its physical origin is the same as for the broad H2*0 component. 
In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear 
that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~ 10 000 AU). The H2O 
line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO 
abundance ratios are measured to be in the range of ~0.1~1, corresponding to H2O abundances of ~10"'-10"* with respect to H2. Approximately 
5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities. 

Key words. Astrochemistry — Stars: formation — ISM: molecules — ISM: jets and outflows — ISM: individual objects: NGC1333 



1. Introduction 

In the deeply embedded phase of low-mass star formation, it is 
often only possible to trace the dynamics of gas in a young stellar 
object (YSO) by analysing resolved emission-line profiles. The 
various dynamical processes include infall from the surround- 
ing envelope towards the central protostar, molecular outflows 
caused by jets ejected from the central object, and strong turbu- 
lence induced within the inner parts of the envelope b y small- 
scale shocks jArce et alj l2007l |j0rgensen et al.1 l2007h . One of 
the goals of the Water In Star-forming regions with Herschel 
(WISH) key programme is to use water as a probe of these pro- 
cesses and determine its abundance in the various components 
as a function of evolution (van Dishoeck et al. in prep.). 

Spectrally resolved observations of the H2O lio-loi line 
at 557 GHz with ODIN and SWAS towards low-mass star- 
forming regions have revealed it to be broad, ~20 km s', in- 
dicative of an origin in shocks (e.g.. lBergin et al.ll2003h . Within 
the large beams (2' and 4'), where both the envelope and the 



entire outflow are present, outflow emission most likely dom- 
inates. Observations and subsequent modelling of the more 
highly excited H2O lines with ISO-LWS were unable to dis- 
tinguish b etween an origin in shocks or an infalling enve- 
lope (e.g.. ICeccarelli et al][T996l iNisini eraDl2002L iMaret et al.l 
2002). Herschel-HIFl has a much higher sensitivity, higher spec- 
tral resolution, and smaller beam than previous space-based 
missions, thus is perfectly suited to a ddressing thi s ques tion. 
Complementary CO data presented by lYildiz et al] (1201 Ol) are 
used to constrain the role of the envelope and determine outflow 
temperatures and densities. 

NGC1333 is a well-studied region of clustered, low-mas s 
star formation at a distance of 235 pc (iHirota et alJ l2008l) . 
In particular, the three deeply embedded, low-mass class 
objects IRAS2A, IRAS4A, and IRAS4B have been ob- 
serve d extensively with ground-based submil limetre telescopes 



(e.g.. lj0rgense n et al.ll2 005l Maret et al. 2005h and interfer ome 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



2OOIL lJ0rgensen et al.ll2007l) . All 
arcmin scales 



ters (e.g., iDi F rancesc o et all 

sources have strong outflows extending over 

(>15 000 AU). Both IRAS4A and 4B consist of multiple pro 



tostars (e.g., [Choi 2005). Because of the similarities between the 
three sources in terms of luminosity (20, 5.8, and 3.8 L©), en- 
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Fig.l. H2O spectra of the t hree NGC1333 so urces. CO 10- 
9 is shown for comparison dYildiz et alJl201Q) : the CO 10-9 
emission in IRAS2A is affected by chopping into outflow ma- 
terial. The top panel shows the decomposition into broad (red), 
medium (blue), and narrow (black) components. The cartoon il- 
lustrates the physical origin of each component. The inset shows 
a zoom on the inverse P Cygni profile in the H2O 2o2-l 11 fine of 
IRAS4A, where the other components have been subtracted; the 
vertical scale ranges from -0.3 to 0.3 K. 



velope mass (1.0, 4.5, and 2.9 Mrr,: |j0rgensen et al.ll2009 t) and 
presumably also age, they provide ideal grounds for comparing 
YSOs in the same region. 

2. Observations and results 

Three sources in NGC133 3, 1RAS2A, 1RAS4A, and 1RAS4B, 
were observed with HIFI (Ide Graauw et alJl2010h on Herschel 
dPilbratt et al.ll2010 !) on March 3-15, 2010 in dual beam switch 
mode in bands 1, 3, 4, and 5 with a nod of 3'. Observations 
detected several transitions of H2O and Hj'^O in the range 
£'u/^B~50-250 K (Table|2]in the online appendix). Diffraction- 
limited beam sizes were in the range 19-40" (4500-9500 AU). 
In general, the calibration is expected to be accurate to ~20% 
and the pointing to ~2". Data were reduced with HIPE 3.0. A 
main-beam efficiency of 0.74 was used throughout. Subsequent 
analysis was performed in CLASS. The rms was in the range 
3-150 mK in 0.5 kms"^ bins. Linear baselines were subtracted 
from all spectra, except around 750 GHz (corresponding to the 
H2O 2ii-2o2 transition) where higher-order polynomials are re- 
quired. A difference in rms was always seen between the H- and 
V-polarizations, with the rms in the H-polarization being lower. 
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Fig. 2. H^**0 lio-loi spectra of the three NGC1333 sources 
along with the CH 536 GHz triplet from the lower sideband (dot- 
ted lines). Spectra are shown for a channel size of 0.25 kms '. 
The spectrum of 1RAS2A has been rebinned to 4 kms ' to il- 
lustrate the detection of a broad feature. The red line shows the 
source velocity at t;LSR~+7.5 kms"'. 



In cases where the difference exceeded 30% and qualitative dif- 
ferences appear in the line profile, the V-polarization was dis- 
carded, otherwise the spectra were averaged. 

All targeted lines of Hj^O were detected and are listed in 
Table [1] and Fig.[T] The lio-loi transition at 557 GHz was not 
observed before the sources moved out of visibility. The H^^O 
lio-loi line was detected in all sources (Fig. |2]i, although the 
detection in IRAS2A was weak (~5cr = 0.13 Kkms '). This 
line is superposed on the ground-state CH triplet at 536 GHz, 
observed in the lower sideband (Fig.|2]l. Neither the Hj'^O ln- 
Ooo nor the 2o2-lii line in IRAS2A is detected down to cr<0.06 
Kkms-i. 

The H2O lines exhibit multiple components: a broad emis- 
sion component (FWHM>20 km s ') sometimes offset from the 
source velocity (yLSR=+7. 2-7.7 kms"'); a medium-broad emis- 
sion component (FWHM~5-10 km s '); and a deep, narrow ab- 
sorption component (FWHM~2 kms"') seen at the source ve- 
locity. The individual components are all reproduced well by 
Gaussian functions. The absorption is only seen in the H2O 1 n- 

000 line and is saturated in IRAS2A and IRAS4A. In IRAS4B, 
the absorption extends below the continuum level, but is not sat- 
urated. Furthermore, the IRAS4A spectrum of the 2o2-lii line 
exhibits an inverse P Cygni profile. The shape of the lines is 
the same within a source; only the relative contribution between 
the broad and medium components changes. For example, in 
IRAS2A the ratio of the peak intensities is ~2, independent of 
the line, whereas in IRAS4A it ranges from 1 to 2. The Hl^^O 
line profiles compare well to the broad component seen in H2O, 
i.e., similar FWHM>20 kms"' and velocity offset. The width 
is much larger than isotopologue emission o f, e.g., C'^O (~l- 2 
km s"') and is centred on the source velocity (lYildiz et al.l2010l) . 
The medium and narrow components are not seen in the Hj^O 

1 iQ-loi spectra down to an rms of 2-3 mK in 0.5 km s"' bins. 

The upper limits to the Hj^O 1 n-Ooo line are invaluable for 
estimating upper limits to the optical depth, t. In the following, 
the limit on t is derived for th e integrated inten sity; in the line 
wings, T is most likely lower dYildiz et al.l 12010 ). In the broad 
component, the limit ranges from 0.4 (IRAS4B) to 2 (IRAS2A), 
whereas it ranges from 1.1 (IRAS4B) to 2.7 (IRAS2A) for the 
medium component of the Hj^O Ih-Oqc line. Performing the 
same analysis to the upper limit on the Hl^^O 2o2-lii line ob- 
served in IRAS2A, infers an upper limit to the optical depth of 
Hj^O 2()2-l 11 of 1.5 for the medium component and 1.9 for the 
broad. Thus it is likely that neither the broad nor the medium 
components are very optically thick. 
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Table 1. H2O and H^**0 emission" in the NGC1333 sources* . 





Transition 


rms^ 




Medium 




Broad 


Source 


^MB 


J Tub dv 


peak 


j Tub dij 






(mK) 


(K) 


(Kkms-') 


(K) 


(Kltms-') 


IRAS2A 


H2O lii-Ooo 


22 


0.22 


2.7 


0.10 


4.1 




2o2-lll 


25 


0.25 


2.5 


0.12 


5.6 




2ii-2o2 


23 


0.16 


1.9 


0.06 


2.7 




3,2-3o3 


79 


0.14 


1.5 


0.06 


2.7 




3,2-221 


52 


0.13 


1.5 


0.08 


3.8 




Hi«OI,„-l„, 


2 




<0.01 


0.01 


0.14 




lii-Ooo 


22 




<0.07 




<0.12 




2(12-1 11 


12 




<0.04 




<0.06 




3l2-3o3 


19 




<0.06 




<0.10 


FWHM''(kms-') 






10.7±0.8 




42±3 


"lsr"* (km s" 


') 




+ 10.7±0.8 




-2.3±3.4 


IRAS4A 


H,0 Iii-Ooo 


23 


0.36 


4.9 


0.35 


15.3 




2o2-lii 


24 


0.34 


3.8 


0.45 


18.0 




2ii-2„2 


23 


0.14 


0.8 


0.43 


14.0 




3l2-221 


100 


0.07 


0.4 


0.32 


13.6 




Hi^Olio-loi 


3 




<0.01 


0.02 


0.43 




lii-Ooo 


23 




<0.06 




<0.()9 


FWHM''(kms-') 






11.1±2.3 




37±4 


"lsr''(1™is- 


1) 






-0.6±0.5 




+8.7±1.0 


IRAS4B 


H20 1ii-O„o 


29 


1.1 


5.1 


0.54 


14.0 




2(12-1 11 


23 


0.63 


3.5 


0.65 


17.6 




2ii-2o2 


17 


0.37 


1.8 


0.40 


10.2 




312-221 


150 


0.32 


1.3 


0.76 


17.5 




H^*0 li„-l„i 


3 




<0.01 


0.02 


0.43 




lii-Oo(i 


16 




<0.03 




<0.04 


FWHM''(kms-') 






4.6±0.5 




24±2 


olsr'' (km s" 


') 






+8.1±0.3 




+8.0±0.5 



Notes. *"' Obtained from Gaussian fits to each component. In the case 
of H2O lii-Ooo, this includes extrapolation over the absorption fea- 
ture. <*' The coordinates used are for IRAS2A: 03:28:55.6; +31:14:37.1, 
IRAS4A: 03:29:10.5; +31:13:30.9, IRAS4B: +03:29:12.0; +31:13:08.1 
(J2000). *'" Measured in 0.5 kms"' bins. '•''' Intensity-weighted aver- 
age of values determined from Gaussian fits of Hj'O emission lines. 
Uncertainties include statistical errors only. 



3. Discussion 

Many physical components in a YSO are directly traced by the 
line profiles presented here, including the infalling envelope and 
shocks along the cavity walls. In the following, each component 
is discussed in detail, and the H2O abundance is estimated in the 
various physical components. 



3. 1 . Line profiles 

The most prominent feature of all the observed line profiles 
is their width. All line wings span a range of velocities of 
~40-70 kms"' at their base. The width alone indicates that the 
bulk of the H2O emission originates in shocks along the cavity 
walls, also called shells, seen traditionally as the standard high- 
velocity component in CO outflow data, but with broader line- 
widths d ue to w ater enhancement at higher velocities (Sect. 13.21 
[Bachiller eFan frggO. Santiago-Garci a et al.ll2009l) . The shocks 
release water from the grains by means of sputtering and in 
high-temperature regions all free oxygen is driven into water. 
The shocked regions may be illuminated by FUV radiation orig- 
inating in the star-disk boundary layer, thus further enhancing 
the water abundance by means of photodesorption. The broad 
emission seen in the H2^0 1 lo-loi line arises in the same shocks 
(see cartoon in Fig.[T]). 

The medium components (FWHM~5-10 kms ') are most 
likely also caused by shocks, although presumably on a smaller 
spatial scale and in denser material than the shocks discussed 
above. For example, the medium component in IRAS2A is seen 
in other grain-product species such as CH3OH (|j0rgensen et alJ 



I2OO5L iMaretet aP I2005L Fig. O, where emission arises from 
a compact region (<1", i.e., <250 AU) centred on the source 
(J0rgensenet al. 2007), and the same is likely true for the 
medium H2O component in that source. In interferometric obser- 
vations of IRAS4A, a small (~few arcsec) blue-shifted outflow 
knot of similar width has been id entified in, e.g., SiO and SO 
(IChoill2005l lj0rgensen et al.ll2007 1. Small-scale structures exist 
in the other sources as well, which may produce the medium 
components. 

The H2O 2o2-lii spectrum of IRAS4A shows an inverse P 
Cygni profile, a clear sign of infall also detect ed in other molecu- 
lar tracers using interferomet er observations dPi Francesco et al] 
I2OOIL Ltergensen et alj 120071) . This infall signature is also ten- 
tatively seen in the In-Ooo line, but here the absorption from 
the outer envelope dominates and little is left of the blue 
emission peak. The signature is not seen in higher-excitation 
lines. The separation of the emission and absorption peaks is 
~0.8 kms~^ wherea s it is ~1.5 kms ' in the observations 
of[D i France sco et alj ( 1200 ll) and larger in the observations by 
J0rgense n et aP (|2007), indicating that the infall observed in 
H2O 2o2-lii takes place over larger spatial scales. 

The passively heated envelope is seen in ground-based ob- 
servations of high-density tracers to produce narrow emission, 
<3 kms"', which may be self-absorbed (Fig.[3]l. For water, this 
type of emission is not seen in any of the sources; the medium 
component is broader by a factor of 2-3 with respect to what 
is expected from the envelope. The absorption seen in all three 
sources is attributed to cold gas in the outer parts of the envelope- 
Using interferometric observations, |j0rgensen & van Dishoeckl 
(120 lOl) detected compact, narrow (~1 kms~') emission in the 
Hj^O 3 13-220 line in IRAS4B possibly originating in the cir- 
cumstellar disk. Scaling the observed e mission to the transi tions 
observed here by assuming 7^1= 170 K (IWatson et alJl2007n . the 
expected emission is typically less than 10% of the rms for any 
given transition. Hence, when extrapolated to the disks surround- 
ing IRAS2A and 4A, the disk contribution to the H2O emission 
probed by HIFI is negligible. The H2O excitation temperature of 
the broad compone nt is 220+30 K, comparable to that found by 
IWatson et alJ (12007 ). but the inferred column dens ity is a factor 
of 100 higher. Thus, the mid-infrared lines seen bv lWatson et al.l 
may come from the same broad outflowing gas found by HIFI, 
provided the mid-infrared lines experience a factor of 100 more 
extinction. 



3.2. Abundances 

3.2.1. Shocks: H2O/CO 

The observed broad compone nts are comp ared directly with 
HIFI observations of CO 10-9 dYildiz et alj|201 0^, because the 
width and position of the lines are similar and they were obtained 
using approximately the same beamsize (22" versus 19"). The 
exception is for IRAS2A, where the blue line wing is not ob- 
served. The advantage is that no detailed models are required to 
account for the H2O/CO abundance, as long as the lines are op- 
tically thin, in particular the emission from the wings. The abun- 
dance ratio is estimated for vario us temperatures by usi ng the 
RADEX escape probability code dvan der Tak et alJl2007) . The 
density is assumed to be 10^ cm~^, appropriate for the large- 
scale core. If the emission is optically thin, the abundance ratio 
scales linearly with density resulting in the same line ratio corre- 
sponding to a higher abundance ratio. There is little variation in 
the predicted ratio for r>150 K, the typical temperature inferred 
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which requires a para-H20 column density of > lO'^ cm"^ if one 
assumes typical values for T and n(H2) of 15 K and 10^ cm~^. 
With a pencil-beam H2 column of ~10^^ cm"^ (J0rgense net alj 
l2002h . the total H2O abundance in the outer envelope is >10 '^. 
For IRAS2A and 4A, the H2O abundance was further con- 
strained using radiative transfer models. The setup is a spheri- 
cal envelope with dens ity and tem perature profiles constrained 
from continuum data (iJ0rgensen et a l. 2009), an infall veloc- 
ity profile V - (2kms"')(r/rin)"'^^, and a Doppler parameter 
b - 0.8 kms '. Line fluxes were computed with the new ra- 
diative transfer code LIME (Brinch & Hogerheijde subm.). The 
models constrained the abundance of water in the outer envelope 
to be -10"*^. Lower values are insufficient to obtain saturated 
absorption in the Ih-Ooo line, and ~10"'^ is the highest abun- 
dance where the resulting narrow emission can be hidden in the 
observed higher-excitation H2O lines. The models predict that 
the H2O emission from the warm inner envelope (r<100 AU) is 
optically thick, hence no constraints can be obtained from the 
H2O spectra on the inner abundance. However, the lack of nar- 
row H^**0 emission infers an upper limit on the H2O abundance 
of ~10"^ (Visser et al. in prep.). 



Fig. 3. Left. Comparison between the medium component in 
IRAS2A and other species observed with ground-based tele- 
scopes. The broad component has been subtracted for easy com- 
parison. The vertical red line indicates the source velocity at 
H-7.7 kms"'. Right. Comparison between H2O 2o2-lii and CO 
6-5 obtained with APEX-CHAMP^, and emission ratios for the 
blue- and red-shifted outflow lobes. 



bv lYildiz etal.l (l2010l) . The line ratios and abundance ratios are 
listed as a function of velocity in Table[3]in the online appendix. 

The abundance ratio increases with increasing velocity from 
H2O/CO of -0.2 near the fine centre to H20/C0>1 in the Une 
wings of all sources for velocity offsets larger than 15 kms"' 
with respect to that of the source (Fig. [3]l. Assuming that the 
CO abundance is 10""*, the H2O abundances are in the range 
of -lO^^-lO""*. Only at high velocities is the temperature high 
enough for oxygen to be driven into water by means of the 
neutral-neutral reaction O + H2 ^ OH + H; OH -H H2 ^ H2O. 
The same result was found in the massive outflow in Orion-KL 
jFranklin et al.ll2008l) . where less than 1 % of the gas in the out- 
flow experiences this high-temperature phase. The fraction of 
gas for which the H2O/CO abundance is >1 is ~5-10% for the 
sources observed here. 

For IRAS2A, a deep spectrum of CO 6-5 obtained with 
CHAMP^ on APEX simultaneously with observations of HDO 
1 1 i-0()o (Liu et al. in prep.) shows the same morphology in terms 
of a broad and medium component (Fig.[3]l. Furthermore, the ve- 
locity offset and FWHM are the same as for H2O suggesting that 
the line profiles are not unique to H2O, although the broad com- 
ponent is far more prominent in H2O. The ratio of peak inten- 
sities for the two components is ~2-3 in H2O versus 10 in CO 
6-5. Analysing the abundance ratio as a function of temperature 
shows that H2O/CO~0.1-l for r>150 K (Table [3]), consistent 
with what is found for CO 10-9. 

3.2.2. Envelope 

The simplest way to constrain the H2O abundance in the outer 
envelope is with calculations using RADEX on the narrow ab- 
sorption in the 1 1 i-Ooo line. The absorption is optically thick - in 
particular for IRAS2A and 4A, where the feature is saturated - 



4. Conclusions 

These observations represent one of the first steps towards under- 
standing the formation and excitation of water in low-mass star- 
forming regions by means of resolved line profiles. The three 
sources have remarkably similar line profiles. Both the H^^O and 
H^**0 lines are very broad, indicating that the bulk of the emis- 
sion originates in shocked gas. The broad emission also high- 
lights that water is a far more reliable dynamical tracer than, e.g., 
CO. Comparing C'**0 to Hj'^O emission and line profiles indi- 
cates that the H2O/CO abundance is high in outflows and low in 
the envelope. Additional modelling of the emission, should be 
able to constrain the total amount of water in the envelope and 
outflowing gas, thus test the high-temperature gas-phase chem- 
istry models for the origin of water. This will be performed for a 
total sample of the 29 low-mass YSOs to be observed within the 
WISH key programme. 
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Table 2. Observed H2O, H'^O and CH transitions" . 



Transition 



(GHz) 



A 



EJkn 
(K) 



(10 



A 

-3 ^-1 



Beam t\ni 
) (") (min.) 



H,0 



lii-Ooo 


1113.34 


269.27 


53.4 


18.42 


19 


43.5 


2o2-lii 


987.93 


303.46 


100.8 


5.84 


22 


23.3 


2ii-2„2 


752.03 


398.64 


136.9 


7.06 


29 


18.4 


3,2-3o3 


1097.37 


273.19 


249.4 


16.48 


20 


32.4 


3,2-221 


1153.13 


259.98 


249.4 


2.63 


19 


13.0 



HfO l,o-loi 


547.68 


547.39 


60.5 


3.59 


39 


64.3 


lii-Ooo'' 


1101.70 


272.12 


52.9 


21.27 


20 


43.5 


2o2-lii 


994.68 


301.40 


100.7 


7.05 


22 


46.7 


3l2-3o3'~ 


1095.16 


273.74 


289.7 


22.12 


20 


32.4 


CH" 3/2,2--l/2,l + 


536.76 


558.52 


25.8 


0.66 


39 




3/2,r-l/2,l + 


536.78 


558.50 


25.8 


0.23 


39 




3/2,r-l/2,0+ 


536.80 


558.48 


25.8 


0.46 


39 





Notes. ^°' From the JPL database of molecular spectroscopy 
( lPickettetaLlll998h . •*' Total on + off integration time. ''' Observed 
in the same setting as the main isotopologue. *''' Observed with H'^0 

lio~loi- 



Table 3. CO 6-5 and CO IO-9/H2O 2o2-lii line ratios in 5 
kms ' intervals and corresponding abundance ratio for r>150 
Kandn=10^ cm"^. 



dULSR 




IRAS2A 




1RAS4A 




IRAS4B 




(kms"') 


CO 6-5/ 


x(H20)/ 


CO 10-9/ 


j:(H,0)/ 


COIO- 


9/ 


c(H20)/ 


COIO- 


9/ 


-((H2O)/ 




H2O202-I11 


A(CO) 


H20 2„2-ln 


x(CO) 


H2O 2o2- 


-111 


x(CO) 


H2O 2o2- 


-111 


.v(C0) 


-20 --15 


5.0 


0.34 






0.8 




1. 11 








-15--10 


3.8 


0.45 






2.0 




0.43 


0.9 




1.00 


-10--5 


4.6 


0.37 






2.8 




0.31 


1.4 




0.64 


-5-0 


9.3 


0.18 






2.4 




0.36 


1.7 




0.50 


0-5 


26.6 


0.06 






2.9 




0.29 


2.3 




0.37 


5-10 


17.0 


0.10 


3.4 


0.26 


3.9 




0.22 


2.9 




0.29 


10-15 


11.1 


0.15 


3.2 


0.27 


3.3 




0.26 


2.1 




0.42 


15-20 


3.5 


0.48 


0.9 


1.00 


2.4 




0.36 


1.6 




0.53 


20-25 


1.4 


1.25 


0.4 


2.00 


I.O 




0.83 


1.1 




0.77 


25-30 










0.8 




1.11 


0.5 




1.67 


30-35 


0.9 


2.00 












0.9 




1.00 



